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Abstract. The resolution associated with the so-called perfect lens of thickness d is 
— 2nd/ ln(|x + 2|/2). Here the susceptibility X i s a Hermitian function in H 2 of the upper half-plane, 
i.e., a H 2 function satisfying x{~ <»>) = xfa). An additional requirement is that the imaginary part 
of x be nonnegative for nonnegative arguments. Given an interval I on the positive half-axis, we 
compute the distance in L°° (I) from a negative constant to this class of functions. This result gives 
a surprisingly simple and explicit formula for the optimal resolution of the perfect lens on a finite 
bandwidth. 
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1. Introduction. The resolution associated with imaging in conventional optics 
is of the order a wavelength. This is a severe limitation in a number of applications 
in nanoscience, e.g., in lithography, microscopy, and spectroscopy. A remedy to this 
impasse — a so-called perfect lens — was proposed by J. B. Pendry [5]. His idea was 
to use metamaterials that allow for tailoring the dielectric permittivity e and the 
magnetic permeability [i by structuring the medium at a length scale much smaller 
than a wavelength pj O [T3] . This may lead to negative refraction and restoration 
of the near- fields. The perfect lens is a negative refraction metamaterial slab of a 
certain thickness d. The metamaterial is considered linear, isotropic, homogeneous, 
and without spatial dispersion. These ideal limits may be approached by appropriate 
metamaterial designs. Note however that certain implementations may inherently 
violate some of these assumptions. For example, in a metal slab there is necessarily 
spatial dispersion (nonlocality) of the dielectric response; this limits the resolution to 
roughly 5 nm [2]. 

As the permittivity and permeability of the perfect lens are negative, the material 
is necessarily dispersive [14]. Thus the perfect lens conditions e = — 1 and /j, = — 1 can 
only be approached at a single frequency. Since most practical applications involve a 
finite bandwidth, this fact limits the performance of the perfect lens. 

In the present work, we will quantify the finite bandwidth behavior. The resolu- 
tion at a single frequency is found to be r — — 2nd/ln(\x + 2|/2), where \ — e — 1 
or x — M — lj depending on the incident polarization. Thus the imaginary parts of 
e and /x and the real parts' deviation from —1 are both crucial for the resolution. 
The formula for r suggests an interesting mathematical problem: Given a physically 
realizable susceptibility and an interval / = [a, b] on the positive angular frequency 
half-axis, find the infimum of \\x + 2||£oo/r\. Quite remarkably, this problem can be 
solved explicitly, and as a result we obtain a simple formula for the optimal resolution 
on a finite bandwidth. 
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Metamaterial 




Fig. 2.1. The perfect lens: A metamaterial slab surrounded by vacuum. The distance from the 
object to the lens plus the distance from the lens to the image equals d. The arrows refer to the case 
when the metamaterial parameters e and fi are close to —I. 



2. Resolution at a single frequency and on a finite bandwidth. Let u) 

denote the (real) angular frequency. The realizability criteria are causality, conjugate 
symmetry, and passivity. They can be expressed as follows: 

X eH 2 , (2.1a) 



= xH, (2.1b) 
Imx(w) > for ui > 0. (2.1c) 

The causality criterion (|2.1ap stating that x belongs to the Hardy space H 2 of the 
upper half-plane, means that the real and imaginary parts of x form a Hilbert trans- 
form pair (Kramers-Kronig relations) jHHjQ The conjugate symmetry (|2.1b[) is a 
result of the fact that the time-domain response function (inverse Fourier transform 
of x) must be real pQ. Condition (|2.1c[) is valid for all passive media, i.e., media in 
thermodynamic equilibrium in the absence of the variable field PQ . 

The resolution as a function of d, e, and [i can be found by solving Maxwell's 
equations OS]. First we assume that the object to be imaged is one-dimensional. The 
slab has orientation orthogonal to the z-axis, the object varies along the x-direction, 
and the polarization of the magnetic field is taken to be along the y-axis, see Fig. 12.11 
The object and the slab are surrounded by vacuum. For each spatial frequency k x 
of the object, we define the transmission coefficient T as the ratio between the plane 
wave amplitudes at the image and the object. By matching tangential electric and 
magnetic fields at both surfaces, we find 



Aek z k' z e l( - k *~ k '^ d 
{ek z + k> z ) 2 e- 2ik '* d - (ek z - k' z f 



T = z (2 2) 



where k z = uj 2 /c 2 — k 2 and k' 2 — e[iui 2 /c 2 — k 2 . The sign of k z must be chosen 
such that Imfc z > 0, while the sign of k' z does not matter in (12. 2|) . We are now 



1 If the medium is conducting at zero frequency, the electric susceptibility e — 1 is singular at 
u) = 0. Then H2.1aH is still valid provided we set x = e — 1 — icr/w, where <r > is the zero frequency 
conductivity JL. In other words, such a medium will have larger loss Ime for the same variation in 
Ree. With no loss of generality we can therefore exclude such media. 
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interested in what happens for large values of k x , corresponding to near-fields that 
decay exponentially in vacuum. We therefore assume that k x S> ui/c. With the 
additional assumptions that |e + 1| <C 1, at most of the order 1, and d at most of 
the order a vacuum wavelength (cod/c < 1), we obtain 

1/T « 1 - [(1 - e 2 ) - (e - n){w/ck x f] 2 exp(2k x d)/l6, (2.3) 

where we have used that 

ek z + k' z = (e 2 k 2 z - kf)/(ek z - k' z ) = [(1 - £ 2 )fc 2 + e(e - ^/cf]/{ek z - k' z ). 

We take the resolution r to be the smallest value of 2ir/k x such that the modulus of the 
second term on the right-hand side of (|2.3p is equal to 1. This definition makes sense no 
matter what this term's phase happens to be, since the exponential factor exjp(2k x d) 
will force \T\ to decrease rapidly when 2ir/k x gets smaller than r. In general, the 
resolution becomes a nontrivial function of e and but if |e — [i\(u> / ck x ) 2 <C 2|e + 1|, 
then 



2nd 



(2.4) 



By the assumption oj/ck x <C 1, the requirement |e — fj,\(w/ck x ) 2 <C 2 je + 1| can be 
rewritten as 



\e- (i\(u/c)d< -2\e+ 1| In 

Thus (|2.4[) is always valid if e = [i. Also, it is valid provided the lens is sufficiently 
thin. Note that in the latter case the resolution is independent of /j,. If we had chosen 
the opposite polarization (i.e., electric field along the y-axis), we would have arrived 
at exactly the same result only with the roles of e and /i interchanged. In other 
words, for a one-dimensional object it is sufficient to have one of the parameters e and 
\x close to —1 with a small imaginary part [5]. For a two-dimensional object, both 
polarizations are necessarily present; thus both e and /i should be close to —1. 

The L°°-norm of the resolution restricted to the interval / measures the poorest 
resolution in the corresponding frequency band. In order to optimize the resolution, 
we should therefore minimize the L°° norm on /. In terms of the electric or magnetic 
susceptibility x = f- lorj; = /'-l ) our task will therefore be to compute the 
following distance: 

inf ||x + 2|| i „ w . 

3. The main result. Before stating our main result, we introduce a few nota- 
tional conventions. 

For < p < oo, the Hardy space H p of the upper half-plane {( = u) + ir\ : 77 > 0} 
consists of those analytic functions / in this domain for which 



/oo 
\f(uj + irf)\ p duj < 00; 
-00 



H°° is the space of bounded analytic functions. A function / in H p has nontangential 
boundary limits at almost every point of the real axis, and the corresponding limit 
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function, also denoted /, is in L p — L P (M.). Indeed, the LP norm of the boundary 
limit function coincides with the H p norm introduced above. Thus we may view H p 
as a subspace of L p . 

As already noted (see (|2.1a[) K we will require the following symmetry condition: 
/(— tj) = f{oS). Functions / satisfying this condition will be referred to as Hermi- 
tian functions. We observe that Hermitian functions have even real parts and odd 
imaginary parts. 

The Hilbert transform of a function u in L p (1 < p < oo) is defined as 

1 f°° «(*) , 
u(u)) — p.v. — / at. 

K J-oo W - t 

It acts boundedly on L p for 1 < p < oo and isometrically on L 2 . If u is a real- valued 
function in L p for 1 < p < oo, then u + iu is in H p , and so the role of the Hilbert 
transform is to link the real and imaginary parts of functions in H p . We will only 
work with Hermitian functions, and we will be interested in computing real parts from 
imaginary parts. For this reason, it will be convenient for us to consider the following 
Hilbert operator: 

Hv(uj) = p.v. - f v(t) ( -J— + - J— ) dt, 

acting on functions in L P (R + ). Provided 1 < p < oo, the function Tlv + iv will then 
be in H p ', with the presumption that v is an odd function. 
For a finite interval I = [a,b] (0 < a < b) set 

_ b 2 - a 2 
~ b 2 + a 2 ' 

Our main theorem now reads as follows. 
Theorem 3.1. 



2A 

inf {\\Hv + k> + 2|| L oo (J) : v > 0, v G L 2 (R+)} 



1 + VI- A 2 ' 

It will become clear that the infimum is not a minimum, but we may extract from 
our proof explicit functions that bring us as close as we wish to the infimum. 

4. Auxiliary results. The main result of [TT] will play a central role in the 
proof of Theorem 13. 11 To state it, we define for each real a the family of functions 

K a (I) = {ue L 2 (R + ) : v{uS) > for w > Q,Hv(w) = a for w G /}. 

(Here and elsewhere we suppress the obvious "almost everywhere" provisions needed 
when considering pointwise restrictions.) We think of functions in K a (I) 7 or more 
generally functions in L 2 (R + ), as the imaginary parts of Hermitian functions, and we 
view them therefore as odd functions on M. 
We also need the function 
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It is taken to be positive for real arguments w > b and is analytic in the slit plane 
C \ ([— b, —a] U [a, b]). For real arguments a < \(\ < b we define a(Q by extending 
it continuously from the upper half-plane. Thus a(C) takes values on the negative 
imaginary half-axis when £ is in (a, b) and on the positive imaginary half-axis when 
—Q is in (a, 6), and otherwise it is real for real arguments. 

Let us also associate with the interval / the following Hilbert operator: 

Hiv{u) = p.v. - f v(t) ( — !— + -—) dt. 

n Jn+\i \t-Lj t + wj 

The main result of [11] was the following parametrization of K a (I). 
Theorem 4.1. A nonnegative function v in L 2 (R + ) is in K a {T) if and only if 
the following three conditions hold: 



v(t)\a(t)\dt < oo (4.1) 

-v 

- { tv(t)a(t)dt = a (4.2) 
Jr+\i 

v(uj)=Hi(av)(uj)/\a(uj)\, ujeI. (4.3) 

The integrability condition (|4.ip is merely a slight growth condition at the end- 
points of /; we may write it more succinctly as 

[v(a -t)+ v(b + t)]^= < oo. 
o yt 

This condition ensures that the integral in (|4. 2[) and the Hilbert transform appearing 
in (|4.3p are both well-defined. 

At first sight, the theorem may not seem to give an explicit parametrization of 
K a (I). However, the Hilbert transform appearing in (|4. 3|) is given by 

If 2t 
Hi (av) (w) = - / »(*>(*)-= %dt, 

7T Jg+y t A - LU 2 

and we observe that the integrand on the right is nonnegative whenever v(t) is non- 
negative. Hence v(w) > for to off I implies v(u) > for u in /. This small miracle 
implies that K a (I) is parameterized by those nonnegative functions v in L 2 (R + \ I) 
for which (14.11) and (14.21) hold and such that 

\Hi (av) (uj)\ 2 \a(uj)\- 2 duj < oo. 

i 

By rephrasing this condition in more explicit terms, we arrive at the following corollary 

Corollary 4.2. A nonnegative function v in L (R + \ /) has an extension to a 
function in some class K a (I) if and only if the following condition holds: 

[u(a - t)u(a -t) + v(b + t)v(b + r)] Ljg!^+l)] dtdr < oo. 

V tr 
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The difference between (I4.1[) and the condition above is the logarithmic factor, which 
means that the condition of the corollary is only a very slight strengthening of (|4.1|) • 
It is clear that for instance boundedness of v near the endpoints of / is more than 
enough. 

We note that the integrand in (|4.2|) is negative to the left of I and positive to the 
right of /. This means that if a is negative, then 

2 f a 

\a\<- tv(t)\a(t)\dt, 



with equality holding if v vanishes to the right of I. It follows that 

Hj M H > - f V (t)a(t)-^^dt > M; (4.4) 

TT Jq t z — UJ Z UJ Z 

we may come as close as we wish to this lower bound by choosing any suitable v 
supported by a small set sufficiently close to 0. 

We remark that the results stated above can be proved by a method similar to 
that given in the next section, the key ingredient being a conformal map sending the 
"two-sided" segment [a 2 ,6 2 ] to the unit circle. 

5. Proof of Theorem l3.1i In what follows, H°° (f2) denotes the space of bounded 
analytic functions on Q equipped with the suprcmum norm, and H P (B>) stands for 
the H p spaces of the open unit disk D. 

Let v be an arbitrary function in L 2 {R + ) such that v > and 

\\Hv + iu||x,w>(j) < ao. 

We will also assumed that v = u|h+w satisfies the condition of Corollary 14.21 We 
know from Theorem 14.11 that the extension of v to a function in some class K a (I) is 
unique, and we may therefore use the notation v for this extension as well. If now tp 
denotes an arbitrary bounded function supported on / with H.<p also bounded, we get 
the inequality 



Setting 



\Hv + iv + 2|| i oo( / ) > inf \\H(f + i<p + iv + 2 + a||zoo(j). 




we may write 

\\Hv + iv + 2|| i oom > inf \\H<p + i(p- -{Hvq - iv ) + 2 + a\\ L <*,m. 

The function H.tp + itp is the boundary limit function of a bounded analytic function 
in the upper half-plane. In fact, since the imaginary part has limit for every point 
in M \ ([—6, —a] U [a, 6]), this function extends by Schwarz reflection to a bounded 



2 This assumption may seem unjustified. However, we may first restrict attention to the smaller 
interval I s = [a + e, b — e] so that v is bounded near the endpoints of I E . Letting e — » 0, we would 
then obtain the same lower bound as we do with our a priori assumption. 
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analytic function / in C* \ ([—b, —a] U [a, b]), where C* = C U {oo}. This function / 
satisfies /(C) = /(C) and /(oo) = 0. The function Hvq + iv extends in the same 
fashion to a function g satisfying g(Q = g(C) and 5(00) = 0. Both these functions are 
in fact analytic in the variable £ = C 2 ; we write 

F(0=f(Vt) and G(£)=giy/t), 

and then it follows that 

w£\\H<p + i<p-l;{Hvo-Wo) + 2 + a\\ L <x,m = inf ^G + 2 + a\\ H <x>r C *\ [a? ,&])■ 

V> 2 F: F(oo)=0 2 

The remaining computation is most easily done if wc first map C* \ [a 2 , b 2 ] conformally 
onto the open unit disk D, say by the map 



b 2 +a 2 \ 2 fb 2 -a 2 ^ 2 



™ = M0 = £nr^U-^r- \/U- 

where the square root is positive for positive arguments. We write T(w) = G(£(w)) 
and obtain 

inf \\H<p + i<p- -^{Hv -iv Q )+2 + a\\ L <*,m = inf \\F + + 2 + a\\ H °°m)- 

tp 2 F: F(0)=0 2 

Since T(0) = F(0) — 0, we have by orthogonality 

\\F+~T+2 + af Hoom > ||r|| 2 ff2(D) + (2 + a ) 2 . 

We may assume a < since for a = 0we may choose T = 0. By (|4.3[) and (|4.4[) . the 
expression on the right-hand side is minimal if 



, \a\y/(b?-u*)(u*-a?) 
vq{u) 5 • 



For I io I = 1, we therefore get 

led ImjD 



Imr(w) 



\- 1 +Rew 
Hence 

2S\a\w 

r H = Y+J^' 

where S < 1 is determined by the equation 

28 



1 + 5 2 

or in other words, 
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It follows that 



1- 



1/2 



F + |r + 2 + a\\ H ~ m > + a) 2 + a 2 ^^) > ^ 

the minimum on the right is obtained when 

a = -2(1 -8 2 ). 

It remains for us to prove the remarkable fact that this lower bound is in fact an 
infimum. To this end, observe that 

2^ + 2(1-^^ = 2^^-. 

1 + ow 1 + ow 

In other words, the minimum is achieved if we choose a = — 2(1 — 5 2 ), tp = 0, and 

2(5(1 -8) 2 w 



T(w) 



Sw 



In view of (|4.3[) and (14. 4[) , we see that we can get as close as we wish to the associated 
minimum by picking a function v > such that i>|r+\j is supported on a small set 
close to 0, 



tv(t)a(t)dt = -2(1 - 6 2 ) : (5.1) 



and 



= -Hi (av) [u)/\a{u)\, w£l (5.2) 

6. Discussion and conclusion. Theorem 13.11 together with (|2.4p gives the op- 
timal resolution as a function of bandwidth: 



2ird 2nd 



111 i+v?-A" ln t 



(6.1) 



where the approximation is valid for A <C 1. In this limit, A w (b — a)/b is the 
relative bandwidth. This optimal resolution is approached when the susceptibility 
x(w) = u(oj) + iv(u>) contains a strong resonance at low frequencies, and a weak 
resonance centered at the relevant bandwidth. As an example, let v(u>) off / be the 
imaginary part of a Lorentzian resonance function, i.e., v(u>) = Im^(w) for ui € R + \I, 
where 

u \ LuJ Q 

l ( UJ ) = — 2^ ■ 

UJq — UJ z — lUJ-f 

Here the resonance frequency oj$ and the bandwidth 7 are positive constants, while 
the strength L is found from the requirement (|5.ip . Computing v(u>) for lu G / using 
(|5.2p . and computing m(u;) = TLv(lu), we obtain the upper plot in Fig. 16.11 Note 
that with this procedure, one may get arbitrarily close to the bound ()6.1|) by choosing 
wo and 7 sufficiently small. Although realizable in principle, it might be difficult to 
fabricate a metamaterial with this near-optimal response. As an alternative, we can 
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Fig. 6.1. The real parts (solid lines) and imaginary parts (dashed lines) of x( w ) + 2 for a 
near-optimal case (upper plot) and a Lorentzian approximation (lower plot). Also shown is the 
associated modulus IxC 1 ^) + 2| (dotted lines). Only the frequency range in the vicinity of I is shown; 
in addition there is a strong resonance atujQ. The Hilbert transforms are computed using fast Fourier 
transforms. The parameters used are a = 2loo, b = 2.1ujo, and ■y = O.Iluq. 

approximate the near-optimal susceptibility by letting x( w ) = u(uj) + iv(ui) be the 
superposition of two Lorentzian functions; one centered at loq and one centered at 
(a + b)/2. We let the first Lorentzian be equal to that in the former example. The 
bandwidth of the second Lorentzian is b — a, and we choose the strength such that 
du((a + b)/2)/du coincides with the corresponding value for the near-optimal case. 
The result is given in the lower plot in Fig. 16.11 

Since the distance from the object to the lens plus the distance from the lens to the 
image equals d, there may be practical reasons for not reducing d. If so, the resolution 
can only be reduced by shrinking the operational bandwidth. Unfortunately, the 
logarithmic dependence of A may require an unpractical, small bandwidth. 

It has been suggested to use a multilayer stack of alternating negative index and 
positive index materials as the lens \±0\. This effectively reduces d in (|6.1jk however, 
then the distance from the object to the lens plus the distance from the lens to the 
image equals the thickness of each layer. In the limit when the layer thicknesses 
approach zero, the resulting effective medium acts as a fiber-optic bundle, but one 
that acts on the near field [lOj . To optimize the resolution and minimize aberrations, 
one again ends up with the problem of minimizing |e + l| and/or If aberrations 

can be tolerated, only Ime and/or Im/x need to be minimized. This can be achieved 
on a finite bandwidth at the expense of some variation of Ree and/or Re/i [Tl I12j. 

We finally note that by simple scaling our result can be used to quantify the 
operational bandwidth of all components with desired permittivity edes and/or per- 
meability /ides less than unity; with a certain tolerance of |edes — e| and/or |/idcs — lA- 
Thus, our result may for instance prove useful for establishing the operational band- 
width of invisibility cloaks [8]. For components with permittivity and permeability 
larger than or equal to unity, there is no bandwidth limitation resulting from (|2.ip . 
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